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We show that the Fermi surface (FS) in the antiferromagnetic phase of BaFe2As2 is composed of 
one hole and two electron pockets, all of which are three dimensional and closed, in sharp contrast 
to the FS observed by angle-resolved photoemission spectroscopy. Considerations on the carrier 
compensation and Sommerfeld coefficient rule out existence of unobserved FS pockets of significant 
sizes. A standard band structure calculation reasonably accounts for the observed FS, despite the 
overestimated ordered moment. The mass enhancement, the ratio of the effective mass to the band 
mass, is 2-3. 

PACS numbers: 71.18.+y, 74.70.Xa 



Precise knowledge of electronic structures near the 
Fermi level Ep is a prerequisite for understand- 
ing the origins of novel superconductivity. In the 
case of the iron-pnictide superconductors pQ, ordinary 
electronic band structure calculations of the parent 
compounds greatly overestimate the antiferromagnetic 
(AFM) ordered moment [2]-[4] . Several moment-reduction 
methods have been proposed, e.g., pnictogen- height 
adjustment [3 ], negative- U [5 j. Moreover, a dynamical 
mean-field theory (DMFT) study suggests that correla- 
tions seriously modify the Fermi surface (FS) in the AFM 
phase [6 . There is a clear need for an experimentally 
determined benchmark FS for judging how successfully 
those theoretical approaches reproduce band structures. 

Recently, it was shown that BaFe2As2 can be de- 
twinned mechanically by compressing or elongating sin- 
gle crystals along a tetragonal [110] axis [7, 8 . This allows 
one to resolve intrinsic in-plane anisotropy in the AFM 
phase. In this Letter, we completely determine the FS 
in the AFM phase of BaFe2As2 via Shubnikov-de Haas 
(SdH) oscillation measurements on detwinned single crys- 
tals. The determined FS is reasonably accounted for 
by our band structure calculation but appreciably differs 
from the FS observed by previous angle-resolved photoe- 
mission spectroscopy (ARPES) studies [9UT3]. We also 
consider implications of the present results for transport 
properties. 

High-quality BaFe2As2 single crystals with the resid- 
ual resistivity ratios of 40-60 were prepared by a self- 
flux method and subsequent anneal [14]. A device sim- 
ilar to one described in Ref. 14 was used to mechani- 
cally detwin them. Standard four-contact resistivity p 
measurements were performed in a dilution refrigerator 
and superconducting magnet. For a sample compressed 
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Fermi 
surface 


Field 
direction 


m* 1 m e 


m* /mband 


r (10- 12 s) 


a hole 


B | 


1 a 


2.8(1) 


3.0 


1.4(1) 




B I 


1 b 


2.7(1) 


3.6 


2.8(4) 




B | 


1 c 


2.1(1) 


2.8 


2.3(2) 


7 electron 


B | 


1 c 


0.9(1) 




2.6(10) 


5 electron 


B | 


1 a 


1.3(1) 


2.0 


1.0(2) 




B | 


1 b 


2.1(1) 


3.0 


2.0(4) 




B | 


1 c 


2.4(3) 


1.7 





(elongated) along a tetragonal [110] axis, the electrical 
current was applied parallel to the compression (elonga- 
tion) direction, resulting in I \\ b (I \\ a) geometry [Fig. 
1(b)]. The electronic band structure of BaFe2As2 was 
calculated within the local spin-density approximation 
(LSDA) using a full-potential LAPW method (TSPACE 
and KANSAI-06). The experimental crystal structure 
at T = 20 K [15] was used. The AFM order of the Fe 
moments was incorporated using the space group Cccm. 
The calculated magnetic moment is 1.6 /i#/Fe. 

Figure 1(a) shows p as a function of temperature T 
near the structural/magnetic phase transition for com- 
pressed, elongated, and freestanding (= twinned) sam- 
ples. The p{T) curves are consistent with previous 
data[3 [8j. Figure 1(c) shows the magnetic-field B de- 
pendence of p in the compressed sample [Fig. 1(b)] for 
B || a, 6, and c at T ~ 0.17 K. After subtracting the 
smooth background, SdH oscillations are clearly visible 
for all field directions [Fig. 1(d)]. The oscillations con- 
tinue to below B — 7 T, much lower than previously 
reported [I6j [17]. Figure 1(e) shows Fourier transforms of 
SdH oscillations in the compressed and elongated sam- 
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FIG. 1: (color online), (a) T-dependence of p near the structural/magnetic phase transition, (b) Schematic diagram of the 
compressed sample, (c) B-dependence of p in the compressed sample for three principal field directions at T ~0.17 K. The 
small magnetoresistance for B \\ b results from the longitudinal configuration I \\ B. (d) SdH oscillations extracted from the 
data in (c) as a function of inverse field 1/B. The curves are vertically shifted for clarity. A third polynomial was fitted to each 
p(B) curve to determine a smoothly varying background pbackground and was subtracted to obtain an oscillatory part p osc - A 
normalized quantity p osc / pbackground is shown, (e) Fourier transforms of the SdH oscillations in the compressed and elongated 
samples for selected field directions. The field angle ca (0 c b) is that between the c axis and B in the ca (cb) plane. The spectra 
are vertically shifted for clarity. The two samples give consistent results (compare the elongated-sample spectra at ca = 49.5° 
and c b = 100.5° with the compressed-sample spectra at ca = 50.5° and B \\ b, respectively). 



pies for selected field directions. The two samples gave 
consistent results. Figure 2(a) summarizes the angular 
dependence of the SdH oscillations in the compressed 
sample. Clear differences are seen between the ca (blue 
squares) and cb (red circles) planes, confirming appropri- 
ate detwinning. 

We find three fundamental frequency branches a, 7, 
and S in Fig. 2(a). The first two branches are consistent 
with (hence named as in) the previous reports p~6j [17], but 
the last branch S has not been observed previously. We 
attribute the previously reported /3 branch to the second 
harmonic of 7 since its angular dependence agrees well 
with the angular dependence of 7. This assignment is 
further supported by the observation that m* associated 
with this frequency is twice that of 7. The a frequency 
shows only small frequency anisotropy, F for B \\ a be- 
ing slightly larger than that for B || 6, but its amplitude 
shows clear anisotropy between the ca and cb planes (note 
the mark sizes in Fig. 2a), which might be explained by 
anisotropic spin-splitting effects [18]. The 7 branch shows 
no detectable in-plane anisotropy. The angular variation 
fits that expected for an FS pocket whose shape is a pro- 
late ellipsoid of revolution, as indicated by the dashed 
line. The angle dependence of the S branch indicates that 
the responsible FS pocket is flattened along the c axis 
and exhibits the most pronounced frequency anisotropy 
between the ca and cb planes. The effective mass m* 



and electron relaxation time r were determined from the 
T- and ^-dependences of the SdH oscillation amplitudes, 
respectively, as usual[19], and are tabulated in Table I. 

Our original band structure calculation indicates that 
bands-68 and -69 cross Ep, giving hole and electron 
sheets of the FS, respectively [Fig. 2(c)]. Each sheet 
consists of large and small pockets. The frequencies orig- 
inating from the large hole and electron pockets are in 
the region of 1 kT, and those from the small pockets are 
in the region of 0.01 kT for B \\ c [pale-colored + and 
x in Fig. 2(b)]. Although the experimental and calcu- 
lated frequencies differ by a factor of approximately two, 
it is clear from the angular dependence that the observed 
a and 6 branches can be assigned to the large hole and 
electron pockets, respectively, of the calculated FS. 

To improve agreement between the experimental and 
calculated a and S frequencies, we shifted the energies of 
bands-68 and -69 by -0.0032 Ry (-44 meV) and +0.0048 
Ry (+65 meV), respectively [Fig. 2(b) thick curves and 
Fig. 2(d)]. Similar procedures have been used in other 
multiband metals very successfully [20, 21 . The agree- 
ment now is satisfactory (our criterion for a is the aver- 
age of the frequencies for the three principal directions), 
though the observed small angular dependence of a is not 
reproduced so well as in the original calculation. The ad- 
justed bands-68 and -69 contain 0.0235 holes and 0.0130 
electrons per primitive cell, respectively [the primitive 
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FIG. 2: (color online) (a) Angle dependence of the experimental SdH frequencies in the compressed sample. The mark sizes 
indicate the oscillation amplitudes. The ca-plane (blue squares) and c6-plane (red circles) data are superimposed. The solid 
lines indicate the fundamental frequency branches. The dotted lines indicating the harmonic and combination frequencies are 
calculated from the solid lines. The dashed line, almost overlapped by a solid line, is a fit to the 7 branch assuming that the 
responsible FS pocket is an ellipsoid of revolution. The fit gives kf = 0.050 A" 1 and k c F = 0.33 A -1 , k F h{c) being the Fermi 
wave number in the ab plane (along the c axis). No frequency other than 3a is observed above F = 1 kT. (b) Comparison 
between the calculated and experimental frequencies. For the calculated frequencies, both original and adjusted ones are shown 
(see text). The solid curves representing the experimental frequencies are the same as those in (a), (c) The FS resulting from 
the original band structure calculation. It consists of a band-68 hole sheet (left) and band-69 electron sheet (right). Since 
our band structure calculation does not include spin-orbit coupling, it predicts chains of alternating hole and electron pockets 
along the grey lines, in accord with Ref. 1231 but they are so thin that they can not be drawn accurately. The directions a, 6, 
and c shown in the figure refer to those of the Fmmm orthorhombic unit cell [15]. (d) The FS in BaFe2As2 determined in the 
present study. The drawings of the a hole and 5 electron pockets are based on the adjusted band structure calculation, and 
the 7 electron pocket is schematically shown, based on the ellipsoid fit in (a). There might be minute FS pockets along the 
grey lines, though our data show no trace of such pockets, (e) Electronic band structure near Ef along TAFZ. The 7 pocket 
is placed at the position marked by the black arrow (see text). 



cell contains two formula units (fu)]. The densities of 
states (DOS) at Ep are 6.64 and 6.35 states/Ry per prim- 
itive cell for the adjusted bands-68 and -69, respectively. 
Comparing the measured effective masses with the band 
masses m^and (Table I), we find average mass enhance- 
ments m* /rriband of 3.1 and 2.3 for the band-68 a and 
band-69 S pockets, respectively; these are small compared 
to the enhancements of 3-20 found in KFe2As2[22 . Us- 
ing these values together with the calculated DOS, we 
can estimate that bands-68 and -69 contribute 1.8 and 
1.2 mJ/K 2 mol-fu, respectively, to the Sommerfeld coeffi- 
cient of the specific heat. 

We now consider the 7 pocket. It is reasonable to as- 
sume that the 7 pocket is located at a position where a 
small pocket appeared in the original calculation. Fig- 
ure 2(e) shows the relevant part of the band structure; 
two likely positions are marked by arrows. The band 



crossings at these points, sometimes referred to as Dirac 
nodes p~2j [23], become anticrossings if spin-orbit coupling, 
which is absent in our calculation, is included and hence 
do not necessarily give FS pockets: Ep may occur in 
a gap. A recent laser ARPES study [24], which is con- 
sidered bulk sensitive because of the low photon energy 
of 7 eV, has however shown that the right (anti) crossing 
marked by the black arrow actually occurs below Ep and 
produces an electron pocket of non- negligible size. We 
identify this pocket with our 7 pocket. The (anti) crossing 
is situated almost at Ep in the original calculation, and 
the renormalized Fermi energy for the 7 pocket is esti- 
mated to be 11 meV from the aforementioned ellipsoid 
fit and measured effective mass. Assuming the mass en- 
hancement to be about three, a band shift of ~30 meV is 
enough to produce the 7 pocket. Because of the symme- 
try, two 7 pockets occur in the BZ [Fig. 2(d)] and enclose 
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0.0116 electrons, resulting, within error, in perfect carrier 
compensation (0.0235 holes vs. 0.0246 electrons). Using 
the measured effective mass for 7, we estimate the contri- 
bution to the Sommerfeld coefficient to be 1.9 mJ/K 2 mol- 
fu. The sum of the contributions from the a, 7, and S 
pockets is 5.0 mJ/K 2 mol-fu, which is in excellent agree- 
ment with a direct measurement on an annealed single 
crystal (5.1mJ/K 2 mol-fu) [25]. The carrier compensation 
and agreement on the Sommerfeld coefficient rule out ex- 
istence of unobserved pockets with a comparable volume 
to the observed ones. Although there might be minute 
pockets with at least an order-of-magnitude smaller vol- 
ume along the grey lines in Fig. 2(d), we have seen no 
trace of such pockets. 

Our band structure calculation accounts well for the 
observed FS if band energy adjustments of at most 65 
meV are allowed, despite the fact that the calculated 
magnetic moment 1.6 \±b overestimates the experimen- 
tal one 0.87 Pb&q\- The magnitudes of the adjustments 
are similar to those necessary in a prototypical multiband 
metal MgE>2 (up to ~ 90meV) [20] . Since our calculation 
is based on the experimental crystal structure and does 
not include ad hoc procedures such as the As-position 
adjustment and negative- £7, those procedures are unjus- 
tified. The observed FS is very different from the FS 
of the DMFT study of Ref. [6[ providing a strong argu- 
ment against the claim that electronic correlations largely 
modify the FS of the AFM phase. 

The agreement between our FS and those observed in 
previous ARPES measurements [9UT3] is limited. For ex- 
ample, the S bright spot, the e petal electron pocket, and 
one of the a and /3 hole pockets reported in Ref. [13] may 
correspond to our 7, J, and a pockets, respectively, but 
the other reported pockets are absent in our FS. Further, 
warped FS cylinders often found in /^-resolved ARPES 
measurements [9j [TTJ [13] are incompatible with our FS 
composed of closed pockets. It should be noted that the 
photoelectron escape depth (~ 5 A) for a typical photon 
energy of 20^40 eV[27] is shorter than the c-axis length 
(13 A), which inevitably limits the k z resolution. 

Our data do not support recent transport studies [28 
attributing observation of the linear-in-5 magnetoresis- 
tance to Dirac fermion transport in the quantum limit. 
Our p(B) curve for B \\ c is clearly superlinear [Fig. 1(c)]. 
The 7 pocket, a supposed Dirac pocket, does not reach 
the quantum limit until B becomes comparable to F 1 = 
90 T. The conductivities of the a and 7 pockets are es- 
timated from m* and r {B \\ c) in Table 1 to be 3.6(2) 
and 5(2) xlO 6 m -1 ^ -1 , respectively. The sum is al- 
ready comparable to the measured conductivity, 11(2) 
xlO 6 m -1 ^ -1 . It is therefore also unlikely that unob- 
served Dirac pockets, if any, dominate the conduction. 

The origin of the in-plane resistivity anisotropy, which 
increases up to pb/p a ~2 with only a few % Co- 
doping [TJ [8], is left elusive. The in-plane mass anisotropy 
is negligible for the a pocket and opposite to the resistiv- 



ity anisotropy for the 5 pocket, i.e., m* along the b axis 
being smaller than m* along the a axis [Table 1, note m* 
for B || a(b) is associated with cyclotron motion in the 
bc(ac) plane]. Intriguingly, this does not seem compatible 
with an optical study [29], which indicates that the Drude 
spectral weight (cj 2 = ne 2 /com in a free-electron model) 
for E || a increases and becomes much larger than that 
for E || b below T/v- It might be better reconciled with 
a more recent study [14] suggesting an isotropic Drude 
component. 

TT thanks Takahiro Shimojima, Tamio Oguchi, and 
Hiroaki Ikeda for helpful discussions. 
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FIG. Angle dependence of the a and 2a oscillation am- 
plitudes. In the cb plane (filled red marks) the ampli- 
tude of a is suppressed near = ±20 and 60°, but no 
such suppression is observed in the ca plane (filled blue 
marks). The observed suppression might be attributed 
to spin-zeros. 1 The measured quantum oscillation is a 
superposition of the oscillations coming from the spin- 
up and the spin-down electrons. They have the same 
frequency but a different phase, resulting in interference 
of the two oscillations. The net oscillation amplitude is 
therefore the single-spin oscillation amplitude multiplied 
by the spin- splitting factor R s = cos(p7rS), where p is a 
harmonic number and S = (l/2)g e jy(m*/m e ). The ef- 
fective ^-factor g e ff has to appropriately incorporate ef- 
fects of many-body renormalization, spin-orbit coupling, 
and antiferromagnetism in the present case. When S = 
1/2, 3/2, the fundamental frequency (p = 1) disap- 
pears (R s = for p = 1) but the second harmonic (p 
=2) remains. If we assume 5 = 1/2 and 3/2 near = 20 
and 60°, respectively, in the cb plane, for example, two 
spin-zeros occur for the second harmonic between these 
angles, which explains the suppressed amplitudes of 2a 
(open pinkish red marks) between the two angles. Since 
m* of the a frequency is almost isotropic in the ab plane 
(Table I of the journal article), the contrasting behavior 
of a between the cb and ca planes has to originate from 
anisotropy of geff] i-e., the observation of spin-zeros only 
in the cb plane seems to indicate that g e ff substantially 
increases as B approaches 6, but not as B approaches 
a. This seems reasonable if we consider the antiferro- 
magnetic order: since the ordered moments are parallel 
to the a axis, the b axis is the direction of the perpen- 
dicular susceptibility of antiferromagnetism, which is not 
suppressed by the ordering, but the susceptibility along 
the a axis is suppressed below Tjy. 
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